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The magnetic properties of the transition metal monoxides MnO and NiO are investigated at equi- 
librium and under pressure via several advanced first-principles methods coupled with Heisenberg 
Hamiltonian MonteCarlo. The comparative first-principles analysis involves two promising beyond- 
local density functionals approaches, namely the hybrid density functional theory and the recently 
developed variational pseudo-self-interaction correction method, implemented with both plane-wave 
and atomic-orbital basis sets. The advanced functionals deliver a very satisfying rendition, curing 
the main drawbacks of the local functionals and improving over many other previous theoretical 
predictions. Furthermore, and most importantly, they convincingly demonstrate a degree of inter- 
nal consistency, despite differences emerging due to methodological details (e.g. plane waves vs. 
atomic orbitals). 
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I. INTRODUCTION 

The relative simplicity of structural and magnetic or- 
dering and the abundance of available experimental and 
theoretical data elect transition metal monoxides (TMO) 
as favorite prototype materials for the ab initio study 
of exchange interactions in Mott-like insulating oxides. 1 
TMO are known to be robust antiferromagnetic (AF) 
Mott insulators with sizable exchange interactions and 
Neel ordering temperatures. The accurate determination 
of magnetic interactions purely by first-principles means 
is a remarkable and as yet unsolved challenged— The dif- 
ficulty stems, on the one hand, from fundamental issues 
in the description of Mott insulators by standard density 
functional theory (DFT) approaches, such as local-spin 
density approximation (LDA) or the generalized gradient 
approximation (GGA). On the other hand, the determi- 
nation of low-energy spin excitations require a meV-scale 
accuracy; however, the error bar due to specific imple- 
mentation and technical differences may easily be larger. 

A large amount of theoretical work for TMO has 
amassed over the years. A number of studies were car- 
ried out in particular for MnO and NiO with a variety of 
advanced methods: the LDA+U^d GGA+U^ifl the 
optimized effective potential (OEP)^ the quasiparticle 
Green function (GW) approach ) 4 ' 11 ! 12 several types of 
self-interaction corrected LDA (SIC-LDA) ,±2r— Hartree- 
Fockr^— and several types of hybrid functionals such 
as B3LYP,22i2a PBEO^i^, Fock-35^ 2 ., and B3PW91.2i 
From a methodological viewpoint, Refs. HHE andflrjlare 



particularly relevant for our present purposes, because of 
the systematic comparison of diverse approaches to com- 
puting magnetic interactions. Other studies^ 7 - - — focused 
in particular on pressure-induced high-to-low spin mag- 
netic collapse observed at very high pressure (~150 GPa 
for MnO) and relevant to Earth-core geophysics. Here we 
will not, however, be concerned with the phenomenology 
of this specific phase transition. 

The present work presents a detailed theoretical anal- 
ysis of MnO and NiO magnetic properties on a wide 
range of lattice parameters (i.e. hydrostatic pressures) 
carried out by an array of both standard and ad- 
vanced first-principles methods. In particular, our the- 
oretical front-liners are two approaches proposed re- 
cently for the description of strongly-correlated systems: 
the Heyd, Scuseria and Ernzerhof (HSE) hybrid func- 
tional approachj22 and the recently developed varia- 
tional pseudo self-interaction correction^ implemented 
in two different methodological frameworks: plane-wave 
basis set plus ultrasoft^ pseudopotentials (PSIC), and 
linear combination of atomic-orbital basis set (ASIC). 
To provide a baseline for their evaluation, we comple- 
ment these methods by their local counterparts imple- 
mented in the same methodological setting, namely LDA 
in plane-waves and ultrasoft pseudopotentials (reference 
for PSIC), local-orbitals and norm-conserving pseudopo- 
tentials (reference for ASIC), and GGA in the Perdew- 
Becke-Ernzerhof^(PBE) version (reference for HSE). 
Performing the same set of calculations in parallel with 
different methods, as well as different implementations 
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of the same method is instrumental to distinguish fun- 
damental and methodological issues, and characterizes 
this work with respect to the many previous theoretical 
studies of TMO. 

MnO and NiO in equilibrium conditions have a high- 
spin magnetic configuration and large (~ 3.5-4 eV) band 
gap. Magnetic moments and exchange interactions de- 
pend crucially on the details of structural and electronic 
properties. The latter are characterized by a complex 
interplay of distinct energy scales: the crystal field split- 
ting, which in rocksalt symmetry separates the on-site 3d 
e g and ti g energies; the charge transfer energy between O 
p and TM d states: the hopping energies between d-d and 
p-d states^ Ref|3 convincingly shows that an empirical 
single-particle potential suitably adjusted to reproduce 
the experimental values of the above mentioned energies 
can deliver highly accurate magnetic interactions (mo- 
ments and spin- wave dispersion). However, obtaining a 
correct balance of all these contributions is difficult even 
for advanced density functional methods, not to mention 
standard LDA or GGA, which fails altogether (to differ- 
ent extents depending on the specific compound). A gen- 
eral analysis of these failures and difficulties of ab initio 
approaches is beyond the scope of this work; we men- 
tion that the thorough analysis carried out in Ref|2| and 
l35l suggests that a single parameter, as adopted by the 
LDA+U, is not sufficient, while global multi-state energy 
corrections could serve this purpose. Both the advanced 
methods (HSE and PSIC/ASIC) adopted in this work, 
while quite different in spirit, act in terms of "global" 
corrections to local density functional energy spectrum, 
i.e. no a priori assumption is made about which particu- 
lar state or band is in need of modification or correction. 

A large portion of the LDA/GGA failure in the de- 
scription of correlated systems can be imputed to the 
presence of self-interaction (SI) - that is, the interaction 
of each single electron with the potential generated by 
itself. One consequence of SI is a severe, artificial up- 
shift in energy of the on-site single-particle energies in the 
density-functional Kohn-Sham equations (compared e.g. 
to measured photoemission energies) . Spatially localized 
states are naturally characterized by a large SI: these 
will include 3d states, but also (albeit to a lesser extent) 
2p states, and of course localized interface, surface, and 
point-defect states. It is sometime argued that this er- 
ror is immaterial because Kohn-Sham eigenvalues are not 
supposed to match single-particle energies, and density- 
functional theories are only valid for ground-state prop- 
erties. However, highest-occupied and lowest-unoccupied 
electronic states are exceptions, since the former must be 
rigorously related to the first ionization energy and thus 
is a ground-state property while the latter may also be 
interpreted within generalized KS schemes [Phys Rev B 
53, 3764, (1996)] as an approximate electron addition en- 
ergy. As a consequence, the SI error may dramatically 
affect the predicted ground state properties, even turn- 
ing the Kohn-Sham band structure from insulating to 
metallic (a frequent occurrence for Mott-like systems in 



LDA/GGA). 

MnO and NiO are both affected by SI, although to dif- 
ferent extents^ severely for MnO, dramatically for NiO. 
Under positive (i.e. compressive) pressure the problem 
will be amplified, as any small band gap which may ex- 
ist at equilibrium will be further reduced up to complete 
closure, and magnetic moments may be disrupted. Thus, 
a reliable description of TMO under pressure necessarily 
requires approaches overcoming the SI problem. Both 
PSIC/ASIC and HSE, although from different starting 
points, work towards the suppression of SI. The former 
explicitly subtracts the SI (written in atomic-like form) 
from the LDA functional; the latter, in a more fundamen- 
tal manner, inserts of a portion of true Fock exchange in 
place of the local exchange functional, whose incomplete 
cancellation with the diagonal Hartree counterpart is the 
source of SI in LDA/GGA functionals<2Z Our results will 
show that, despite the different conceptual origin, the two 
approaches deliver a systematically consistent description 
of MnO and NiO, in fact with spectacular quantitative 
agreement in several instances. 

The paper is organized as following: Section |TT] de- 
scribes briefly the methodologies employed; in Sec- 
tion IIII1 the model used to calculated the exchange- 
interaction parameters is discussed. Section [IV] illus- 
trates our results at equilibrium (IIV Ap and under pres- 
sure (|TVB"1) . In Sections HVBl] and UVU we discuss the 
exchange interactions and critical temperatures, respec- 
tively. Section IVl offers some concluding remarks. 

II. METHODS 

Our first-principles results are obtained by using three 
different codes: PWSIG^, SIESTA 3 ^ an d VASPiS. With 
the PWSIC code, which uses plane-waves basis set 
and ultrasoft pseudopotentials we carry out calcula- 
tions within LDA (hereafter LDA-PW) and PSIC. The 
SIESTA 3 ^ package, implemented over a local (atomic) 
orbital basis set and norm-conserving pseudopotentials, 
was employed for LDA (LDA-LO below) and ASIC cal- 
culations (note that the ASIC is available only in a devel- 
oper version of the code). The variational PSIC/ASIC 31 
approach evolves from the non-variational precursors 
PSIG2&ii and ASIC42 The new method describes elec- 
tronic properties as accurately as its precursors, and vari- 
ationality allows its application to structural optimiza- 
tion and frozen-phonon calculations via quantum forces. 
In the following we briefly sketches the main features of 
the method, remanding the reader to Ref. [3l| for the de- 
tailed formulation for both plane waves (PSIC) and local 
orbital (ASIC) basis set. 

The PSIC starts from the following energy functional 

ijva 

where here i, j are collective quantum numbers (/,-, m^, 
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lj=h, raj) of a minimal basis (i.e. one for each angu- 
lar moment I) of atomic wavef unctions, and a and v are 
the spin and the atomic site index, respectively. Thus, 
to the usual LDA energy a SI contribution is subtracted 
out. This is written as single-particle SI energies (e SI ) 
scaled by the effective orbital occupations (p). For an 
extended system, whose eigenfunctions are Bloch states 
(Vrik)' ^ ne orbital occupations are self-consistently calcu- 
lated as projection of the Bloch states onto the atomic 
orbitals 6 



(2) 



In addition, the "effective" SI energies are defined as 
4- = E /5k <Ckl7<,*>C« (7,>Kk) (3) 

nk 

where ^f i v is the projection function of the SI potential 
for the i th atomic orbital centered on atom v. In the 
radial approximation one has 

7i,,wii,f( r ) = VHXc[Pu,h(r); 1] <f>k,rrH,v( r ) > ( 4 ) 

with Cij being normalization coefficients 

c-l 



= / dT(f>i umi (r) V H xc[Pu,u(' r ) : > l ] < t > k,m j {*) -(5) 



It is now easy to see that the functional derivative of 
this PSIC energy functional leads to Kohn-Sham single 
particle Equations similar to those adopted in the original 
scheme described in Ref. [ill ]. 

The VASP— code, emloying the projected augmented 
wave (PAW) approach, is used for PBE^ 3 . and HSE^ii cal- 
culations. In the HSE method the short-range (sr) part of 
the exchange interaction (X) is constructed by a proper 
mixing of exact non-local Hartree-Fock exchange and ap- 
proximated semi-local PBE exchange. The remaining 
contributions to the exchange-correlation energy, namely 
the long-range (lr) exchange interaction and the elec- 
tronic correlation (C), is treated at PBE level only, re- 
sulting in the following expression: 
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(6) 

The partitioning between sr and lr interactions is 
achieved by properly decomposing the Coulomb kernel 
(1/r, with r = |r — r'|) through the characteristic param- 
eter n, which controls the range separation between the 
short (S) and long (L) range part 

L W + M , ) = ^W + *1. (7 ) 

r r r 

We have used here /i = 0.20 A -1 , in accordance to the 
HSE06 parameterization^ 3 - and corresponding to the dis- 
tance 2/yit at which the short-range interactions become 



negligible. For /i— 0, HSE06 reduces to the parent un- 
screened hybrid functional PBEO 4 ^. In terms of the one- 
electron Bloch states ^kn(r) and the corresponding occu- 
pancies /kn the sr Hartee-Fock exchange energy E 1 ^'^'^ 
can be written as 



E 



HF,sr,fj, 
X 



E 2w k/kn X Wqf< 
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(8) 



(9) 



x ^ n (r)^ qm (r)^ m (r')0kn(r') , (10) 



from which one can then derive the corresponding non- 
local sr Hartee-Fock exchange potential™. In the last 
few years the application of the HSE06 to a wide class 
of solid state systems, including archetypical benchmark 
examples 4 ^, transition metal oxide a 21 ' 47 i 48 , dilute mag- 
netic semiconductors 4 ^^, multifferoics^i and magnetic 
perovskites^, has demonstrated, that HSE06 delivers a 
substantially improved description of ground state struc- 
tural, electronic, magnetic and vibrational properties 
with respect to the standard PBE and PBE+U. How- 
ever, the use of hybrid functionals for metals in less sat- 
isfactory, because of the general significant overestimated 
bandwidth 4 ^. 

As far as the technical aspects are concerned, PWSIC 
calculations have been carried out by using 16-atom face 
centered cubic (FCC) supercells [i.e. 8 formula units 
(f.u.)], cut-off energies of 40 Ry, reciprocal space integra- 
tion over 6x6x6 and 10x10x10 special k-point grids for 
self-consistency and density of states calculations, respec- 
tively. VASP calculations have been performed using a 4 
f.u. unit cell, an energy cut-off of 25 Ry, a 4x8x4 k-point 
mesh and a standard HSE mixing parameter a=0.25. 
The SIESTA calculation have been performed using a 4 
f.u. cell with a 6x11x6 k-point mesh, a real space mesh 
cut-off of 800 Ry and a double-^ polorized basis. Pres- 
sures have been evaluated using the Birch-Murnaghan 
equation of stated 

Montecarlo simulations of the classical 2-parameter 
Heisenberg model have been carried out for a spin lat- 
tice system of size L=12 (i.e. N = L 3 total lattice sites). 
We determined ground state magnetic ordering and crit- 
ical temperature by simulated annealing for each pair of 
ab-initio-calculated Ji and J2 parameters characterizing 
the magnetic structure (see next Section), at each lattice 
constant and for each method. In order to tests finite-size 
effects on the results some annealing with L=20 was also 
performed. The annealing was done over 30 temperature 
points, starting from high temperature (roughly twice the 
critical temperature as estimated by few trial runs) down 
to T=0, with 10 6 sweeps at each temperature. The an- 
nealing protocol is the usual Metropolis algorithm based 
on single spin update. 
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III. MAGNETIC STRUCTURES AND THE 
HEISENBERG MODEL 



TMO have a rock-salt structure (see FigO}, so each 
TM has 12 nearest-neighbor (NN) and 6 next-nearest 
neighbors (NNN). The NNN are connected through oxy- 
gen bridges, and their interaction J 2 is dominated by 
superexchange. On the other hand, NN interact via a 
typically smaller exchange coupling Ji whose sign may 
depend on the specific TMO; Ji involves direct TM- 
TM exchange (giving a robust AF contribution) and a 
90°-oriented TM-O-TM superexchange (expected to be 
weakly FM). The observed ground state magnetic phase 
is antiferromagnetic (111) A-type, labeled AF2 hereafter. 
It can be seen as a stacking of (111) planes of like spin 
alternating along the [111] direction, as illustrated in Fig. 
[TJ In AF2 each TM has 6 spin-paired intra-(lll)-plane 
NN and 6 spin-antipaired inter- (lll)-plane NN; on the 
other hand, all 6 NNN bonds are inter-planar and an- 
tipaired. Thus, this configuration maximizes the energy 
gain associated to the NNN antiparallel spin alignment. 
As for beyond-NNN magnetic interactions, there is am- 
ple experimental^ and theoretical^ evidence that they 
can be safely discarded (e.g. according to inelastic neu- 
tron scatterings in NiO they are two order of magnitude 
smaller than the dominant J2. We explicitly checked this 
with the SIESTA code.). 

In order to evaluate Ji and J2 we need to consider 
at least two competing high-symmetry magnetic phases 
beside the observed AF2. Natural choices are the fer- 
romagnetic (FM) order and the AF (110) A-type order 
with (110) spin-paired planes compensated along [110] 
(labeled AFi). AFi can also be seen as made of FM (001) 
planes alternating along [001] (see Fig. [IJ. The AFi 
phase has all the 6 NNN spin-paired, 4 of the NN spin- 
paired, and 8 NN spin-antipaired. On the other hand, 
it is interesting to note that the G-type AF order (also 
depicted in Fig. [1} is strongly disfavored by frustration, 
since in FCC symmetry there is no way to arrange the 12 
NN interactions in antiparallel fashion without conflict. 

To extract Ji and J 2 we fit our calculated total energies 
to a standard 2-parameter classical Heisenberg Hamilto- 
nian of the form: 



H 



Ji ^ el ■ e~j -J2 ^] el ■ e~j 



(11) 



where and {{i,j)) indicate summation over NN 

and NNN, respectively, and el is the spin direction unit 
vector. Energies (per f.u.) are then expressed as: 



Efm = E — 6 Ji — 3 J2 
Eafi = E + 2Ji ~ 3J2 
Eaf 2 = E + 3 J2 




0=r 




FM 




AFI: A-type (001)/(110) 



T ho 




AF2: A-type (111) 



G-type 



FIG. 1: (Color on-line) Magnetic phases used for the Heisen- 
berg model fit: FM; AF2 that is built of alternating (111) 
planes of like spins [highlighted by thick (red) lines]; AFi, 
made of alternating (011) (or equivalently (001)) planes of 
like spins, delimited by (red) thick lines. The highly frus- 
trated G-type AF phase is also shown for comparison. Filled 
light (green, up spin) and black (down spin) circles indicate 
TM atoms, small empty circles are for oxygens. For the G- 
type phase (blue) filled circles marked with a question mark 
indicate TM atoms with frustrated spin direction. 



This is solved to give: 



Ji — $(Eaf! - E F m) 

Ji = ^(4:Eaf 2 — 3Eaf! — Efm) 



(13) 



(12) 



With this choice of the Hamiltonian, negative and 
positive J values correspond to energy gain for spin- 
antiparallel and spin-parallel orientations, respectively. 

Finally, we mention that several other anisotropic 
terms may in principle contribute to the Heisenberg 
Hamiltonian, related to short-range dipolar interactions 
favoring a preferential spin direction parallel to (111) 
planes, and to rhombohedral distortions of the AF 2 phase 
consisting on a (111) inter-planar contraction and slight 
change of the perfect 90° angle of the rock-salt cell, which 
causes a symmetry breaking of Ji in two and J]~ 
values~A However, all these effects are quantified to be 
order-of-magnitude smaller than the dominant exchange- 
interaction energies (e.g. for NiO jf -J^ ~ 0.03 meV ac- 
cording to neutron data's), thus the Heisenberg Hamil- 
tonian written in Eq llll can be considered fully sufficient 
for our present purposes. 
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TABLE I: Equilibrium lattice constants ao (in A) and bulk 
moduli Bo (in GPa) calculated in this work with various meth- 
ods, compared with experimental values. All values refer to 
the stable AF2 magnetic ordering. 

LDA-PW PSIC LDA-LO ASIC PBE HSE Expt. 
MnO 



ao 


4.38 


4.35 


4.35 4.35 


4.47 


4.41 


4.43" 


Bo 


158 


194 


177 199 


145 


170 


151 6 ,162 c 








NiO 








ao 


4.15 


4.09 


4.09 4.06 


4.19 


4.18 


4.17 d 


Bo 


234 


269 


270 287 


183 


202 


180-220 e 


a): 


Ref.5,5; b): 


Rcf.56. 


c): Ref.57 








d): 


Ref.58; e): 


Ref.59 









IV. RESULTS 
A. Equilibrium structures 

We have calculated total energies and pressures of 
MnO and NiO as a function of lattice parameter for the 
3 magnetic phases FM, AF 2 , and AFi. Values of the 
equilibrium lattice parameter and bulk modulus for the 
stable AF2 phase are reported in Table HI in comparison 
with the experimental values. 



Results are quite satisfactory overall: each method pre- 
dicts an equilibrium lattice constant in good agreement 
(within 1-2%) with experiment for the AF2 phase. It is 
well known that structural properties calculated by LDA 
or GGA can be good, or even excellent, although the elec- 
tronic properties are poori 45 i 47 The results for MnO and 
NiO are a case in point, as both LDA-PW and PBE stay 
within 1% from the experiment (the former in defect, the 
latter in excess), while the LDA-LO result is about 2% 
below the experimental value. As for beyond-local func- 
tionals, HSE sllightly underestimates PBE results and as 
a consequence is quite close to experiment, while PSIC 
and ASIC underestimate their respective local-functional 
(LDA-PW and LDA-LO) references by ~1% on average 
(a tendency also found^L in other classes of oxides such 
as titanates and manganites). 

In Figs. [5] and [3] the calculated pressures for, respec- 
tively, MnO and NiO are reported as a function of lattice 
parameters for the three magnetic phases. The behavior 
in the region around the equilibrium lattice constant is 
expressed by the calculated bulk modulus (Bo) in Tab|l] 
The advanced functionals coherently give an increase of 
Bo by ~20-30 GPa (~15%) with respect to their respec- 
tive local functionals. This increase is a consequence of 
the enhanced 3d state localization and concomitant in- 
crease in Coulomb repulsion under compression which 
is expected from beyond-local approaches. Concerning 
the agreement with experiment, both PBE and HSE give 
values within the reported experimental uncertainty. In 



contrast, LDA-PW gives Bo at the higher end of the ex- 
perimental error bar, thus that the 15% further increase 
caused by PSIC pushes the value of B ^40-50 GPa above 
the experiment. Hence the discrepancy should be seen 
as due to the LDA-PW performance (and to the charac- 
teristics of the used pseudopotentials), rather than as a 
failure of the PSIC method in itself. This trend can also 
be observed for the local orbital basis, where the LDA 
results for Bo are already slightly above the the experi- 
mental values, and the addition of the ASIC corrections 
simply pushes the bulk modulus further away. 

A very interesting feature which emerges consistently 
from all the methods is the quite similar pressure de- 
pendence for different magnetic orderings, especially ev- 
ident for NiO. This looks surprising at first glance, espe- 
cially at strong compressive pressure, where (as discussed 
later on) changes in the magnetic ordering are related to 
metal-insulating transitions and to radical changes in the 
electronic properties. The explanation is that e g elec- 
trons (strongly hybridized with O p states) govern the 
electronic and magnetic properties, but have only a mi- 
nor effect on the response to applied hydrostatic pressure. 
In MnO, where ti g states are also magnetically active, 
pressure is slightly more sensitive to the specific mag- 
netic ordering (in FigJ5]AF2 tends to differ from AFi and 
FM, which almost overlap each other). It should also be 
noted that in MnO for very contracted lattice constants 
the advanced functionals give pressures a factor of 1.5-2 
larger than those of the corresponding local functionals, 
depending on the method and specific magnetic phase. 
At variance, for NiO advanced and local functionals give 
pressures in the same range. This reflects the larger ef- 
fect of the advanced functionals on the half-filled ti g shell 
of MnO, which is pushed down in energy and increase 
its spatial localization and its Coulomb repulsion under 
compression, than on the filled ti g shell of NiO. 

B. Magnetic properties upon applied pressure 

1. Magnetic phase diagram under pressure 

Figs. U and [S] summarize our findings concerning phase 
stability and magnetic moments for MnO and NiO. Each 
panel reports results obtained by a given energy func- 
tional for relative magnetic energies (with respect to the 
most stable magnetic ordering) and their correspond- 
ing magnetic moments, as a function of lattice constant. 
Column-wise, panels are ordered according to the code 
used: LSDA-PW and PSIC results (left, PWSIC); PBE 
and HSE results (center, VASP); LSDA-LO and ASIC 
results (right, SIESTA). 

We start our analysis from MnO results given by LDA- 
PW, LDA-LO, and PBE (top panels of Figg}. We 
can capture immediately the substantial similarity of the 
three approaches: at large volumes all have a region 
where AF2 is the most stable phase and insulating, a 
metallic AFi region ^50-100 meV above AF2, and a 
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FIG. 2: (Color on-line) Calculated pressure for MnO in the three considered magnetic phases (see text). Each panel reports 
results obtained by a different energy functional or methodology considered in this work. 




FIG. 3: (Color on-line) Calculated pressure for NiO in the three considered magnetic phases (see text). Each panel reports 
results obtained by a different energy functional or methodology considered in this work. 



metallic FM region ^100-150 meV above AF 2 . Going 
from large- to small-lattice region, all three methods re- 
ports a phase transition (indicated by the vertical dashed 
lines) through which the AF2 phase yields to a metallic 
FM region as ground state, which wins over higher-energy 
AFi and (further above) AF2 metallic phases. The criti- 
cal value of lattice parameter corresponding to the tran- 
sition differs in the three approaches 4.3 A in LDA-PW, 
4.15 A according to PBE, and 3.9 A for LDA-LO. Look- 
ing at the corresponding magnetic moments, LDA-PW 
and PBE show a similarly gradual moment decay when 
going from 4.7 A up to a threshold of 4.0 A (LDA-PW) 
and 4.1 A (PBE), corresponding, in the pressure scale, 
to~40 GPa and 80 GPa, respectively. Below this thresh- 
old, magnetic moments fall quite abruptly down to ~ 2 
(is at 3.8 A for the stable FM metallic phase (notice 
that this unphysical collapse described by LDA-PW and 
PBE is not related to the true moment collapse^ in MnO 
at much higher pressures). At variance, moments cal- 
culated by LDA-LO decay gradually through the whole 
lattice constant range, without collapsing as in LDA-PW 
and PBE. This difference reflects the persistence of the 
AF2 stability in a larger parameter range described by 
LDA-LO. 

Now we analyze the results obtained within beyond- 
local functionals. Overall, the picture is radically differ- 
ent: for the whole range of lattice parameters, the AF2 
insulating phase is robustly the ground state, and the 



spurious phase transition discussed above is absent. Fur- 
thermore, all three approaches coherently report a sta- 
bility enhancement (i.e. a roughly linear energy gain) 
of AF2 for decreasing lattice constant. This effect is 
indeed expected as a consequence of the increased Mn 
d-0 p covalency and the related strengthening of AF su- 
perexchange coupling. The AF2 maximum stability is 
reached at ^3.9 A according to PSIC (corresponding to 
an applied compression of ~130 GPa), at ^4.0 A for 
ASIC (P~130 GPa) and HSE (P~90 GPa). The peak of 
AF2 energy gain with respect to the equilibrium struc- 
ture is nearly 100%, from ^50 meV/f.u. to more than 
100 meV/f.u. according to HSE and PSIC. Above AF 2 
all our advanced functionals favor the AFi phase, that, 
at variance with the always insulating AF2 phase, un- 
dergoes a metal- insulating transition at 3.97 A (PSIC), 
3.95 A (ASIC) or 4.05 A (HSE). Above AFi resides a 
FM region, again separated in a large- volume insulating 
and small-volume metallic sides, with an insulating-metal 
transition threshold of 4.07 A for PSIC, and 4.15 A for 
both ASIC and HSE. This consistency is also reflected in 
similar values of magnetic moments through the whole 
lattice constant range: all methods describe a moderate 
decline from ~ 4.7 fis at 4.7 A to ~ 4.0-4.2 fig (depend- 
ing on the specific magnetic phase) at a=3.8 A. 

Interestingly, the remarkably coherent picture deliv- 
ered by our advanced functionals for MnO is not limited 
to the predict the same ground-state, but also involves 
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Lattice constant (A) Lattice constant (A) Lattice constant (A) 

FIG. 4: (Color on-line) Total energies (relative to the ground state) and magnetic moments of FM, AFi and AF2 phases of 
MnO as a function of lattice parameter, calculated by LDA-PW and PSIC (left panels), PBE and HSE (center panels), LDA-LO 
and ASIC (right panels). The insulating (I) or metallic (M) character of each magnetic phase is also indicated. Vertical dashed 
lines indicate phase transitions. 



ordering and energy differences among the three mag- 
netic phases. This is instrumental to coherently describe 
finitc-tcmpcrature properties as well, as shown in Section 
IIVB3I 

Now we move to the analysis of NiO results, summa- 
rized in Fig. [5j starting again from the phase stability 
diagram drawn by the three local functionals (upper pan- 
els). At variance with what seen for MnO, now all the 
methods give the insulating AF 2 phase as stable at any 
lattice constant; however the competition with the other 
two orderings is described differently: according to LDA- 
PW, moving from large to small lattice constants there 
is first a tiny region where the AF2 stability increases, 
reaches maximum at 4.4 A (thus much above the equi- 
librium value) and then start falling linearly all the way 
down to 3.8 A . Furthermore, above AF 2 the LDA-PW 
predicts a coexistence of degenerate AFi and FM metallic 
phases. This scenario can be rationalized looking at the 
magnetic moments for AFi and FM calculated in LDA- 
PW: starting from the large lattice constant value of ~ 1 
fi b , the magnetic moment falls rapidly and vanishes alto- 
gether just at 4.4 A (i.e. still above the equilibrium value 
of 4.35 A). Below this threshold the LDA-PW describes a 
metallic, Pauli paramagnetic region. On the other hand, 
the expected Mott-insulating behavior is only maintained 
in the AF 2 phase. This feature represents a major short- 
coming which seriously hamper the NiO description by 
LDA-PW. 

PBE shows a similar, although slightly less dramatic 



failure, since the moment collapse starts to show up for 
smaller lattice constant values (4.0 A for FM and 4.1 
A for AFi, thus definitely below the equilibrium 4.19 
A ), and the magnetic moment is severely reduced to 
about 0.5 fiB, without vanishing completely. Quite dis- 
tinct is the behavior of LSDA-LO, which shows two well 
energy-separated phases, AFi above AF 2 , and then FM 
further above, both undergoing metal-insulating transi- 
tions when moving from smaller to larger lattice parame- 
ters. This peculiarity is clearly reflected in the calculated 
magnetic moment: there is no collapse at least until 4.0 
A, then a rapid downfall starts, but the moments remain 
quite sizable up to 4.0 A , at variance with what seen for 
PBE and LSDA-PW, and consistently with what is found 
for MnO as well. This signals the tendency of LDA-LO 
to conserve sizable Mott gap and magnetic moments in a 
range of lattice parameters where plane- wave based im- 
plementations describe the system as already collapsing 
to metallic paramagnets (as explicitly shown later on in 
the DOS analysis). 

Moving to the analysis of beyond-LDA functionals, all 
of them consistently predict the large enhancement of 
AF 2 stability upon lattice constant decrease in a wide 
range around the equilibrium value. The artificial mo- 
ment collapse described by the local functionals is absent, 
and all the magnetic phases remain insulating through 
the whole lattice constant range. All methods find mag- 
netic moments of about 1.7-1.8 fiB at large lattice con- 
stants, and a very moderate decrease to ~1.5 [1b at the 
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FIG. 5: (Color on-line) Total energies (relative to the ground state) and magnetic moments of FM, AFi and AF2 phases of NiO 
as a function of lattice parameter, as calculated for the different methods used in this work: LDA-PW and PSIC calculations 
(left), PBE and HSE calculations (center), LDA-LO and ASIC calculations (right). PM indicates Pauli Paramagnetic ordering, 
I and M insulating and metallic character, respectively. Vertical dashed lines indicate phase transitions. 



smallest lattice constant considered (3.8 A). Particularly 
striking is the agreement between HSE and PSIC, both 
describing a tiny FM region intermediate between AF2 
and AFi, and a parabolic rise of the AF2 energy gain 
from -50 meV at 4.7 A up to ~250 meV at 3.8 A. On 
the other hand, ASIC gives a different phase ordering 
(AF 2 -AFi-FM), with a smaller increase of the AF 2 rel- 
ative stability from ^0 meV at 4.7 A to ^120 meV at 
a=3.8 A. 

In order to clarify the difference in the magnetic mo- 
ments under pressure obtained by the different methods, 
we have calculated the orbital-resolved density of states 
(DOS) for NiO at two representative lattice constants, 
a— 4.5 A and a=4.0 A , corresponding to situations of 
tensile and compressive strain. The results are shown in 

Fig. m 

Starting from the local-functional results at 4.0 A, 
we notice a radical difference between the two LDA 
implementations. LDA-LO predicts a metallic ground 
state for both the FM and AFi phases, with sizeable e g 
spin-splitting and substantial magnetic moments on Ni, 
whereas the AF2 phase is found insulating, with a ^1 
eV energy gap. On the other hand according to LDA- 
PW FM and AFi phases are Pauli-paramagnetic, with 
perfectly compensated spin-densities, while in the AF2 
phase a barely visible gap opens up within the e g man- 
ifold. PBE results stand somewhat in the middle of the 
previous two descriptions: for FM and AFi some Pauli 
magnetization shows up in the upper e g manifold, while 



for AF2 the magnetic moment is formed, and a Mott gap 
is about to open. Even at a=4.5 A some differences 
among the local functionals are detectable: according to 
LDA-LO a minimal gap is opened for each phase, while 
for LDA-PW only the AF2 is insulating. PBE provides 
well formed magnetic moments in each phase but only 
AF2 is clearly insulating. 

Going to the results for beyond-local function- 
als, the expected picture of wide-gap intermedi- 
ate charge-transfer/Mott insulator described by the 
experiments^— is restored. Now the DOS is that of a 
robust insulator under both lattice expansion and com- 
pression, with a valence band top populated by a mixture 
of O p and Ni d states, and the conduction bottom with 
a majority of e g and a minor fraction of O p states. The 
energy gap for the AF 2 phase (~3 eV for ASIC and ^3.5 
eV for both PSIC and HSE) is in good agreement with 
the experimental value, and even FM and AFi phases 
show sizeable gaps of about 1-2 eV. Overall, the similar- 
ity of spectral redistributions for e g , t% g , and O p states 
(especially for HSE and PSIC) is remarkable. 

In summary, for both MnO and NiO beyond-local func- 
tionals deliver a very coherent description of relative 
phase stabilities in the whole examined range of lattice 
parameters, and predict a clear enhancement of the AF2 
phase relative stability (not described by local function- 
als) within a wide lattice parameter interval, which sug- 
gests the possibility of an enhancement of the magnetic 
ordering Neel temperature (T^v) upon applying compres- 
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FIG. 6: (Color on-line) Orbital-resolved DOS for NiO calculated at two lattice constant values (4.0 A, and 4.5 A) with all our 
employed functional: LDA-LO (top left), ASIC (top right), LDA-PW (middle left), PSIC (middle right), PBE (bottom left) 
and HSE (bottom right). Only the DOS for the relevant orbitals are shown: O p (filled gray curves) and Ni d, separated in 
t2 g (filled red curves) and e 9 (solid blue lines) contributions. Positive and negative curves represent majority and minority 
contributions. 



sive stress. Before exploring the validity of this expecta- 
tion we will first discuss the evolution of the magnetic 
coupling constants upon compression. 



2. Exchange interactions under pressure 

A few qualitative considerations on magnetic interac- 
tions can help to correctly interpret our results. Our 



calculations find the T=0 magnetic ground state to be 
the observed AF2 for both MnO and NiO; however, the 
detailed magnetic interactions suggest two different sce- 
narios. Magnetic coupling between Mn 2+ 3d 5 ions is me- 
diated by half- filled orbitals (thus Ji mainly by ti g -ti g 
and J 2 by e g -e g couplings), which are both robustly AF 
oriented according to superexchange theory.— >§2 Thus we 
expect Ji and J2 to be both sizable and negative (i.e. AF 
in our present convention). For Ni 2+ d s ions, on the other 
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hand, only e g -e g coupling is magnetically active. Hence 
we expect a large and negative J2 due to the dominance 
of covalent superexchange, very small and positive Ji due 
to superexchange-mediated 90°-oriented e g -(0 p x ,Py)-£g 
orbital coupling, and therefore huge J2/J1 values. 

These expectations are largely confirmed by our re- 
sults: in MnO (Fig|TJ) Ji and J2 roughly track each other 
as function of the lattice parameter irrespective of the 
calculation method. However, it is of the utmost im- 
portance to observe the dramatic difference between the 
description of local and advanced functionals: looking 
at LDA-PW results [panel a)] Ji and J2 are moderately 
negative at expanded lattice, then upon lattice shrink- 
ing they both change sign and grow up to a maximum 
value at ~ 4 A, and finally fall back down as the lattice 
shrinks further. The PSIC almost completely reverses 
this behavior: Ji and J2 nearly vanish at 4.7 A (sig- 
naling a shorter interaction range with respect to LDA) , 
and then grow steadily (in absolute value) on the nega- 
tive side as the lattice squeezes up. Curiously, LDA-PW 
and PSIC curves intersect each other at ~ 4.4 A, but this 
agreement near the equilibrium lattice is just a fortuitous 
crossing of two otherwise radically different behaviours. 
Notice finally that while in PSIC J2/Ji>l at any lattice 
parameter, in LDA-PW the exchange interaction ratio 
fluctuates as function of the lattice constant. 

The considerations exposed for LDA-PW and PSIC 
can be identically repeated for PBE and HSE, respec- 
tively (the similarity of curves is apparent comparing 
panel a) and b) in FigJT]). The behavior of LDA-LO 
[panel c)] is more peculiar: as previously noticed, the 
local-orbital implementation seems to cure, though only 
partially and perhaps accidentally, already at the LDA 
level some of the deficiencies of the local functionals re- 
lated to the self-interaction problem. It turns out indeed 
that the exchange interactions determined by LDA-LO 
show characteristics similar to those depicted by the non- 
local functionals: at large lattice parameter both Ji and 
J2 are vanishing and negative, then they correctly grow 
negative upon lattice shortening up to 4.2 A(thus well 
below the equilibrium structure), but then they revert 
slope and turn positive below 4.0 A. The ASIC, on the 
other hand, restores through the whole lattice range the 
correct description already seen for PSIC and HSE. 

Now we move to examine NiO (Fig[S]). As expected, 
the relative weight of Ji and J2 is very different: all 
functionals (both local and beyond- local) describe Ji as 
very small and positive at any lattice value (black sym- 
bols). On the other hand, the large and negative J2 is 
again differently described by the two sets of function- 
als. For LDA-PW/PBE and PSIC/HSE (FigJHJ we can 
repeat most of the considerations made for J2 in MnO: 
the behavior with lattice parameter is roughly inverted, 
although the two curves crosses at different values of 
the lattice constant (4.2 A for LDA/PSIC, 3.9 A for 
PBE/HSE). LDA-LO is again a case in its own, as it ap- 
pears to provide results not far from the ASIC ones in 
the lattice region above 4.2 A. ASIC qualitative behavior 
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FIG. 7: (Color on-line) Exchange-interaction parameters Ji 
(black filled circles) and J2 (red filled squares) as a func- 
tion of lattice parameter calculated for MnO with various 
approaches. From top to bottom: a) with PWSIC code by 
LDA-PW and PSIC functionals; b) with VASP code using 
PBE and HSE functionals; c) with SIESTA code using LDA- 
LO and ASIC functionals. Dashed lines refer to local den- 
sity functional (LDA-PW, PBE, LDA-LO) calculations, solid 
lines represent beyond-local density functionals calculations 
(PSIC, HSE, and ASIC). Black and red crosses show exper- 
imental values for Ji and J2, respectively, reported for MnO 
in Table [UJ 

for both Ji and J2 is in line with that of PSIC and HSE, 
although its J2 appears amplified by nearly a factor 2. In 
consideration of what we have seen so far, the reason of 
this difference seems to reside more in the methodology 
(i.e. technical implementation: plane wave versus atomic 
orbitals) rather than in the employed functional. 

From a phenomenological viewpoint, it is important 
to note the appreciable growth of the exchange inter- 
actions for decreasing lattice parameter in a wide re- 
gion extending well below the equilibrium value (i.e. for 
large applied compression), coherently described by all 
our beyond-local functionals. In the next section we will 
illustrate the remarkable reverberations of this behavior 
on the magnetic ordering temperature. For now, we re- 
mark that this is the expected behavior of the so called 
covalent exchange, i.e. the shorter the TM-0 distance, 
the stronger the energetic advantage for the O 2p lig- 
and states to overlap with the adjacent TM 3c? states 
with unlike spins. This advantage reaches a maximum at 
a certain compression, after which the J's start falling. 
This is the point when the pressure is so strong that mi- 
nority DOS begins to be appreciably populated and in 
turn magnetic moments start falling. 

Overall, the comparative analysis of the three beyond- 
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lattice parameter (A) 

FIG. 8: (Color on-line) Exchange-interaction parameters Ji 
(black filled circles) and J2 (red filled squares) as a function of 
lattice parameter calculated for NiO with various approaches. 
From top to bottom: a) with PWSIC code by LDA-PW and 
PSIC functionals; b) with VASP code using PBE and HSE 
functional; c) with SIESTA code using LDA-LO and ASIC 
functionals. Dashed lines refer to local density functionals cal- 
culations (LDA-PW, PBE, LDA-LO), whereas solid lines rep- 
resents beyond-local density functionals calculations (PSIC, 
HSE, and ASIC). Black and red crosses show experimental 
values for Ji and J2, respectively, reported for NiO in Table 

im 

local functionals is very satisfying, as they furnish a quali- 
tatively and quantitatively coherent and physically sound 
description of exchange interactions under pressure for 
MnO and NiO. There is a large body of data in litera- 
ture, both theoretical and experimental, with which to 
compare our data, at least at equilibrium. In Table ILT1 we 
report our values for the theoretical equilibrium structure 
in comparison with the experimental values and other 
theoretical predictions obtained by several density func- 
tional based approaches (see Section I; we neglect tight- 
binding or shell-model results, which rely on experimen- 
tal fitting). 

The J's reported in the Table are determined by finite 
energy differences, unless otherwise specified. In some 
cases, the magnetic force theorem (MFT) 64 based on the 
exchange-correlation density-functional gauge invariance 
under infinitesimal spin rotation a 2 ' 65 was used. In Refill, 
exchange interactions and the whole spin-wave spectrum 
was determined from the poles of spin susceptibilityj 66 ! 67 
As for experiments, results from both inelastic neutron 
scattering (INS)^^ an d thermodynamic data (TD)£&Z2 
are reported. 

Comparison of J values given by different approaches 
and considerations on the level of agreement with the 



available experimental data must be taken very carefully. 
Differences of the order of a few meV may derive from 
technical implementation aspects rather than from the 
basic theory itself, as testified by the sizable difference 
of results obtained by the same theory in different cal- 
culations (e.g. our LDA-LO and LDA-PW). Concerning 
MnO, the comparison is further complicated by the close 
similarity of the two J values. It is expected that local 
functionals should overestimate the J's, due to the un- 
derestimation of intra-atomic exchange splitting and the 
related overestimation of p-d hybridization. This is in- 
deed verified with LDA-LO and PBE. Notice, however, 
the relevant exception of LDA-PW (overall, the least ac- 
curate of the three local functionals) which deliver J's in 
excellent agreement with the experiment: this fortuitous 
agreement was previously explained as consequence of 
the unphysical AF2 to FM transition occurring in LDA- 
PW near the equilibrium structure. On the other hand, 
our beyond-local functionals perform quite satisfactorily, 
apparently ranking among the closest to experiment both 
in absolute terms and concerning the J2/J1 ratio (1.16 for 
INS data, 1.1 for HSE, 1.4 and 1.5 according to ASIC and 
PSIC). 

For NiO the analysis is simpler, as J 1 is very small and 
a qualitative comparison can be done on the base of J2 
only. We have already commented that the LDA-PW is 
grossly inadequate, thus it can be left aside. As for LDA- 
LO and PBE, they deliver sizably overestimated J2. As 
for beyond-local functionals, we have previously seen that 
that J's calculated by PSIC and HSE almost overlap each 
other throughout the lattice parameter range. The agree- 
ment with values drawn from neutron experiments^ is 
indeed quite satisfying. Nevertheless, the slight PSIC un- 
derestimation of the equilibrium lattice parameter (4.09 
A against the near experimental-matched 4.18 A of HSE) 
reverberates in a ^15% overestimation of J2. On the 
other hand, ASIC delivers J's that barely differ from 
the corresponding LDA-LO values, and are amplified by 
nearly a factor 2 with respect to PSIC and HSE. Looking 
at previous literature, we found a substantial agreement 
of PSIC and HSE with GGA+U calculations of Refi 
(here the J's are also calculated as a function of lattice 
parameter) and with other types of hybrid functionals 22 
as well. On the other hand, both unrestricted HF 22 , full 
SIC in LMTO approach (SIC-LMTO)^ and the local 
SIC (LSIC) (a KKR-based implementation of the self- 
interaction correction methocU^— ) tend to an excessive 
electronic localization, which thus turns into a slight un- 
derestimation of the exchange interactions. 



3. Critical transition temperatures under pressure 

Fig. reports critical temperatures for MnO and 
NiO calculated with the Heisenberg Hamiltonian given 
in Eqlll] and solved through classical MonteCarlo (MC) 
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TABLE II: Exchange interaction parameters for MnO and 
NiO (in meV) calculated in this work, compared to experi- 
mental and theoretical values from previous works. 

MnO NiO - 
Ji J2 Ji J2 

Experiment 



INS a -4.8 -5.6 INS 6 1.4 -19.0 

TD C -5.4 -5.9 TD d -1.4 -17.3 





This work: local functionals 






LDA-PW 


-2.7 -6.3 LDA-PW 


-0.5 


-14.7 


LDA-LO 


-10.6 -13.6 LDA-LO 


5.4 


-43.5 


PBE 


-9.5 -14.9 PBE 


1.2 


-44.5 




This work: advanced functionals 






PSIC 


-5.0 -7.6 PSIC 


3.3 


-24.7 


HSE 


-7.0 -7.8 HSE 


2.3 


-21.0 


ASIC 


-8.0 -11.3 ASIC 


5.2 


-45.0 



Previous calculations 



LSDA e (MFT) 


-13.2 


-23.5 








LDA+IT (MFT) 


-5.0 


-13.2 


GGA+U' 


1.7 


-19.1 


OEP e (MFT) 


-5.7 


-11.0 


SIC-LMTO m 


1.8 


-11.0 


PBE+U 7 


-4.4 


-2.3 


Fock-35™ 


1.9 


-18.7 


PBE0 / 


-6.2 


-7.4 


B3LYP" 


2.4 


-26.7 


HF / 


-1.5 


-2.32 


UHF™ 


0.8 


-4.6 


QPGW 9 


-2.8 


-4.7 


QPGW fl 


-0.8 


-14.7 


LSIC' 1 


1.4 


-3.3 


LSIC h 


2.8 


-13.9 


LSIC' 1 (MFT) 


-1.8 


-4.0 


LSIC h (MFT) 


0.3 


-13.8 


B3LYP ,; 


-5.3 


-11.0 









a): Ref.68, b): Ref.54, c): Ref.69 d): Ref.70 
e): Ref.2, f): Ref.21, g): Ref.4, h): Ref.Lo 
i): Ref.23, 1): Ref._9, m): RefT5, h): Ref.22 



simulated annealing technique^ Values calculated at 
equilibrium and at experimental volume are reported in 
Tab. IIII1 in comparison with the experiment and other 
theoretical predictions. 

We start discussing the case of MnO as described by 
our local density functionals. As expected, MC calcu- 
lations describe fictitious phase transitions (consequence 
of the spurious magnetic moment collapse previously dis- 
cussed) from AF2 to FM metallic magnetic phase (high- 
lighted by the dashed areas) while moving from large 
to small lattice parameters. The onset of this transi- 
tion depends on the method: it is especially harmful in 
LDA-PW as it occurs near (just below) the theoretical 
lattice value (4.35 A); much less damaging in LDA-LO, 
where it happens at a much larger lattice contraction 
(4 A); PBE is halfway the two previous cases. In all 
cases this phase transition dramatically alters the criti- 
cal temperature behavior, which is expected to grow (at 
least in some interval around equilibrium) as the lattice 
parameter contracts. Furthermore in all cases the pre- 
dicted Tn (see Tab. IIIII) remains quite distant from the 
experimental Tn=118 K (horizontal dashed line in the 



Figure). Notice again that while LDA-LO and PBE over- 
estimates T^r at the equilibrium structure, the spurious 
phase transition in LDA-PW cause T^r to be smaller than 
the experimental values, and paradoxically not too far 
(at experimental lattice constant) from the experiment. 
This is another manifestation of the fortuitous agreement 
already pointed out in the illustration of the exchange in- 
teractions. 

Conversely, the advanced functionals deliver for MnO a 
nicely consistent picture, with T n growing steadily from 
very large lattice constants (Tjv —>•()) up to 3.9-4.0 A, and 
peaking at -410 K (PSIC and HSE) or 480 K (ASIC). 
The corresponding pressure is around 120±20 GPa de- 
pending on the approach (see Figj^]). There is a minor 
offset between the three methods, due to the slight in- 
crease in J's moving from PSIC to HSE to ASIC. How- 
ever, if calculated at their respective equilibrium values, 
both PSIC and HSE predict a Tjv which is almost spot- 
on to the experimental value. 

Moving to the analysis of NiO, the three local func- 
tionals behave quite differently from each other, although 
PBE and LDA-LO fortuitously cross nearby the equilib- 
rium value. So, it is all the more remarkable that their 
beyond-local counterparts are capable to rebuild a very 
consistent picture, with Tn linearly growing along with 
the decrease of lattice parameter. The near-overlap of 
HSE and PSIC is especially striking, and already noticed 
for the J's. The ASIC curve appears translated upward 
by nearly a factor 2, coherently with the overestimate of 
the dominant J2 value. Notice that for NiO there is no 
lattice parameter turning point in the considered range, 
thus Tj\r keeps growing up to 3.8 A, corresponding to a 
pressure of about 100-120 GPa (see FigJ3]). The agree- 
ment with the experiment for NiO is much less outstand- 
ing than for MnO, as both HSE and PSIC remains below 
the experimental K by —20% (at the equilib- 

rium structure the PSIC value is actually not too far from 
the experiment thanks to its underestimated lattice pa- 
rameter), while ASIC is close to the much overestimated 
LDA-LO value. 

The disagreement between ab initio theoretical esti- 
mates and the experimental value of Tjv for NiO is typ- 
ical and well documented in literature (see e.g. Refflq 
and references therein) and a full discussion on the sub- 
ject is beyond our present scope. We only remark that 
the mismatch with HSE and PSIC is somewhat puzzling, 
in consideration of the excellent agreement of the calcu- 
lated J2 with the experimental value. In fact, even using 
experimental sets of J's, the MC-calculated T^v would 
change only marginally our theoretical value: as a use- 
ful reference we also report in the Table the Tn of MnO 
and NiO calculated by MC using sets of experimental J's 
reported in Table HU It turns out that these T n under- 
estimate by ~ 30% the directly measured TV (they are 
even lower than those obtained with our calculated J's, 
since the slight overestimation of our calculated J's helps 
in shifting up the predicted Tjv). The discrepancy be- 
tween experimental J's and Tjv somewhat points out to 
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lattice parameter (A) lattice parameter (A) 



cluster exact diagonalization (CED). The latter seems to 
restore an outstanding agreement with the experimental 
Tat for NiO£2. 



TABLE III: Critical temperatures (K) for MnO and NiO 
calculated in this work at equilibrium and experimental (in 
brackets) lattice constant, compared with experimental (neu- 
tron diffraction, ND) and theoretical values from previous 
works. As a reference we also report MC-calculated values 
obtained by using the experimental J's (see text for discus- 
sion) . 

MnO NiO 
Experiment 

ND a ' 6 118 523 

Expt. J's 85 c 340 d 

Expt. J's 90^ 300^ 

This work: local functionals 

LDA-PW 55 (71) 272 (280) 

LDA-LO 220 (176) 967 (851) 

PBE 249 (257) 824 (829) 

This work: advanced functionals 

PSIC 116 (89) 458 (387) 

HSE 125 (116) 393 (400) 

ASIC 182 (137) 1048 (850) 

Previous works 



FIG. 9: (Color on-line) Critical temperatures as a function 
of lattice constant for MnO (top panel) and NiO (bottom) 
calculated by simulated-annealing MonteCarlo simulation of 
the Heisenberg Hamiltonian in Ea llll Left and right panels 
separate local and advanced functionals used to determine 
each set of (Ji, J2). The shaded areas in the top-left panel 
indicate FM metallic regions; apart from that, each curve 
separates low-T insulating AF2 from high-T insulating PM 
regions (PM stands for Pauli Paramagnetic ordering). The 
dashed horizontal lines indicate experimental Tat values. 



LSIC 9 (DLM) 126 336 

LSIC 9 (RPA) 87 448 

LSIC 9 (MC) 90 458 

LDA' 1 (MC) 423 965 

LDA+U^MC) 240 603 

CED' 1 (MC) 172 519 



a): Ref.73, b): Ref.74, c) Ref.68, d) Ref.54, e) Ref.69. 
f) Ref.70, g) Refl6, h) Ref3. 



possible inadequacies of the employed Heisenberg model, 
possibly due to further terms (not included in EqlTTj) 
which might be important at the relatively high ordering 
temperature of NiO. 

In Tab lHIl we compare our results with some previ- 
ous theoretical predictions (we omit the many mean-field 
approximations which are known^ to grossly overesti- 
mate the critical temperature). Ref. UM proposes Tjv 
obtained by disordered local moments (DLM), random- 
phase approximations (RPA) and MC, based on the 
MFT-calculated J's given in Table HI The DLM values 
are closest to our HSE and PSIC estimate, while RPA 
and MC values for NiO are larger than our values despite 
smaller J's, due to the debatable inclusion in Ref. 16 of the 
quantum rescaling factor (S+1)/Sj^ Refill also proposes 
MFT-calculated J's, derived from LDA, LDA+U, and 
LDA plus dynamical mean field approach solved through 



V. SUMMARY AND CONCLUSIONS 

It is fair to affirm that the overall account of the 
structural, electronic, and magnetic properties of MnO 
and NiO provided by the advanced functionals is over- 
all quite satisfying, internally consistent, and in good 
agreement with experiments. In particular, HSE shows a 
remarkable quantitative agreement with experiments on 
most examined properties; the PSIC, perhaps surpris- 
ingly when considering the substantially different con- 
ception at the basis of their theoretical constructions, is 
quite comparable with HSE results, and in some cases 
in spectacular quantitative agreement (e.g. the NiO ex- 
change interactions vs. lattice constant). An important 
persistent shortcoming of PSIC or ASIC, however, is the 
prediction of structure: the predicted lattice constant is 
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below experiment by ~l-2%. This tendency to deliver 
smaller-than-optimal structural parameters was also en- 
countered in other situations^! and it is probably not an 
isolated occurrence, but rather a characteristics of the 
PSIC/ASIC method. In perspective, we expect that this 
drawback could be overcome by adopting the GGA (e.g. 
PBE) instead of the LDA as reference functional upon 
which to build the PSIC projector (see Refs l36ll4ll for de- 
tails). This would probably lead, as in the case of HSE, 
to a moderate volume reduction compared to the slightly 
overestimated GGA volume, hence probably to an end 
product much nearer to experiments. 

Some additional considerations concern the compari- 
son of different implementations of the same theory, i.e. 
LDA-PWvs. LDA-LO and PSIC vs. ASIC. Interestingly, 
some visible quantitative differences appear already at 
level of local functionals. Hence, these are not induced 
by the ASIC/PSIC formulation, which in fact tends to 
re-equalize the respective descriptions. The amplitude of 
the PSIC vs. ASIC discrepancy is only minor for MnO, 
and a bit more pronounced for NiO, especially concerning 
the amplitude of J2. This is not particularly surprising, 
since different choices of local orbital basis set and/or 
pseudopotentials can easily alter meV-scale quantities. 
In fact we can find in the literature a wide spread of 
values for exchange interactions predicted by different 
implementations of the same theory (even at the LDA 
level). Not by chance the different basis set has more 
impact on NiO exchange interactions than in the MnO 
ones, since the e g charge spreads more in space than the 
highly localized t2 S one, and as such it is clearly more 
sensitive to the incomplete description of the intersticial 
space due to local orbital basis. As a trade off, this larger 
error bar is compensated by the ASIC capability to tackle 
large-size systems (up to a few thousands atoms) which 
in fact are impracticable for the, in principle, more ac- 
curate plane-wave approaches such as PSIC and HSE. 
A valuable testing ground of this capability may be the 
field of oxide interfaces, where a plethora of new excit- 
ing phenomena have been recently discovered. It will be 
interesting to evaluate whether PSIC/ASIC or HSE are 
capable to describe e.g. composite systems such as an in- 
terface between a metal and insulator, where the degree 
of charge localization at the two sides is rather different. 

Finally, our evaluation of the exchange-interaction pa- 



rameters and of the Neel temperatures requires a men- 
tion. While the calculated J's of both MnO and NiO 
are found in satisfying agreement with the experiments, 
only for the first an equally satisfying T^ is predicted by 
the MC-solved Heisenberg Hamiltonian. This discrep- 
ancy can stimulate debate and more work devoted to 
investigate the possible inadequacies of the Heisenberg 
Hamiltonian at high temperature, an aspect that has not 
been sufficiently stressed in previous literature. 

In conclusion, we have presented a systematic anal- 
ysis of the structural and magnetic properties of MnO 
and NiO under applied pressure, including ground-state 
and finite temperature properties, by using a range of 
standard and advanced first-principle approaches. The 
advanced techniques (HSE, PSIC, ASIC) describe very 
consistenly the behavior of the exchange interactions in 
a wide range of lattice constant values around the equi- 
librium structure, showing an overall agreement with ex- 
periments. This places such methods among the most 
accurate now available to the first-principles community. 
Our results establish a benchmark of accuracy for in- 
novative techniques aimed at the determination of the 
magnetic properties of magnetic oxides. 
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